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AFEUi-NASA  CoUaboration 


ri 


Provide 

economical,  accurate  sensitivities 
for  multidisciplinary  design 
and  analysis 


Computational  AeroSciences 
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SAMS  Overview 
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Cost  &  Capability 


i 


Concept  Technology  System  Development  Production  &  Operations 

Refinement  Development  &  Demonstration  Deployment  Support 


Ignorance  &  Uncertainty 


Cost  of  change 


%  LCC  locked  in  %  LCC  incurred 

•  Knowledge  is  nnost  limited  when  majority  of  cost  is  locked  in 

•  Limited  knowledge  inevitably  leads  to  mistakes 


•  Changes  become  more  expensive  over  time 

•  Correcting/accepting  mistakes  means  cost  growth  and/or  lost  capability 


Better  early 
modeling 


Increased 

knowledge 


Fewer 

mistakes 


Successful 

design 


informed 

decisions 


DISTRIBUTION  STATEMENT  A:  Approved  for  public  release.  Distribution  is  unlimited. 


4 


Technology  Impact 


'A 


'a 


Candidate 

Technologies 


Engineering 

Capability 


Mission 

Capability 


Aeroelastic  wing 
Flutter  suppression 
Third  stream  engine 
Advanced  materials 
Thrust  vectoring 
Innovative  control  effectors 
Directed  energy 

Conformal  load  bearing  antennas 


Weight 

Drag 

Roll  performance 
Yaw  performance 
RCS 

Fuel  efficiency 
Sensor  efficiency 


•  Life  cycle  cost 

•  Survivability 

•  Maintainability 

•  Kills  per  $ 

•  IDs  per  $ 


^(Mission  Capability) 

Use  accurate  physical  models  to  obtain - “Z — : - : - : - 

^  ^  a(Technology) 
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Contribution 


•  CFD-based  analytic  aeroelastic  sensitivities 

•  Steady  and  transient  solutions 

•  Nonlinear  structural  dynamics 

•  Demonstrated  capability 

•  Enhanced  FUN3D 
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Components 


LEGEND 

I  SAMS  component 
External  component 
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Schedule 


YEAR  1  YEAR  2  YEAR  3 


Initiation 

Modularization 
Integration 
Mode-based  analysis 
FEA-based  analysis 
Finalization 
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Current  and  Potential  Customers 


Government  institutions 


Academia 


0 


PrattaWfcitnBy 


Aerospace  companies 


Other  industries 
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Hermes 
Integration  Framework 
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Layered  Networking 


Components  in  a  multidisciplinary  simulation  each  have  their  own 
communication  requirements,  driving  a  layered  networking  model 

MSTC  Engineering  «<or  similar»> 

•  MADO  process  model 

•  Geographically  distributed  components 

Hermes  «<built  for  SAMS»> 

•  Fast  coupling  for  unsteady  analysis 

•  Locally  distributed  components 

MPI 

•  Very  fast  coupling  within  a  component 

•  Execution  on  a  cluster/supercomputer 


Notional  diagram 
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Hermes  Framework 


Characteristics 

•  Remote  procedure  calls 

•  Cross-language 

-  C++ 

-  Fortran  2003 

-  Python 

•  Built  on  0MQ 

-  Client-server  communication 

-  Message  based 

-  Several  transport  protocols 

-  Several  messaging  patterns 


Remote  Procedure  Call 


CLIENT 

SERVER 

appantrt 

AddIi  cation  - —  Procedures 

II 

♦  t 

- IT — — 

call 

\  t 

Client  Stub 

Server  Stub 

\  t 

\  t 

Client 

Server 

Runtime  Library 

Runtime  Library 

\  f 

\  t 

Transport 

Transport 

t 
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Process 


Object 

files 


Executable 


Hermes 

library 


LEGEND 

I  Written  by  a  developer 
I  Hermes  components 
Other  tools 
Generated  files 
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Interface  Definition 


structure  point  { 
real64  x, 
real64  y 

} 

exception  oops  { 
int32  code 

} 

interface  geometry  { 

point  centroid(vector<point>  ps)  throws  (oops), 
real64  distance(point  p,  point  q) 


} 


Defines  a  server’s  remotely  callable  procedures 
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Server  and  Client  Code 


Server 

from  geometry  import  Point,  Geometry 
import  math 
import  zmq 

ciass  Server(Geometry.Server): 
def  distance(seif,  p,  q): 
dx  =  q.x  -  p.x 
dy  =  q.y-p.y 

return  math.sqrt(dx  *  dx  +  dy  *  dy) 
def  centroid(seif,  ps): 

X  =  sum([p.x  for  p  in  ps])  /  ien(ps) 
y  =  sum([p.y  for  p  in  ps])  /  ien(ps) 
return  Point(x,  y) 

addr  =  ‘tcp://*:49152’ 

ctx  =  zmq.ContextO 

srv  =  Server(ctx,  addr,  zmq. REP) 

srv.serveO 

srv.cioseO 

ctx.destroyO 


Client 

from  geometry  import  Point,  Geometry 
import  zmq 


addr  =  ‘tcp://iocalhost:49152’ 
ctx  =  zmq.ContextO 

cii  =  Geometry.Ciient(ctx,  addr,  zmq.REQ) 

p  =  Point(1.2, 3.4) 
q  =  Point(5.6,  7.8) 
c  =  cli.centroid([p,  q]) 
d  =  cii.distance(p,  q) 

cii.cioseQ 

Ctx.destroyO 


Developer  can  focus  on  the  business  logic 
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Future  Directions 

•  Expand  type  support  in  Fortran 

•  Develop  better  documentation 

•  Optimize  for  performance 

•  Support  more  programming  languages 

•  Asynchronous  ’’futures”  model 
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MAST  &  YAMSS 
Structural  Solvers 
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MAST 


Multidiscipliiiciry  Design,  Adaptation,  &  Sensitivity  Toolkit 


•  Elasticity,  heat  transfer,  &  compressible  flow 

•  Adjoint  solver  for  sensitivity  analysis 

•  High-order  finite  elements 

•  Feature-based  and  goal-oriented  adaptation 

•  MPI-  and  thread-based  parallelism 
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Structures 


•  Buckling,  modal,  steady,  &  transient  problems 

•  Linear  and  nonlinear  elements 

•  Geometric  and  material  nonlinearities 

•  Can  solve  linear  analysis  about  nonlinear  point 
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Example  Application 


transient  nonlinear 
thermoelasticity 


panel  design  for  transonic  flutter 
about  nonlinear  equilibrium 
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Analyze  Across  Disciplines 


/ 


Transient  Heat  Transfer 


Sensitivity  Analysis 

Jvy/  db 


Flutter 


flutter  velocity,  Vj- 
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Solution  Process 


•  MAST  is  a  C++  programming  library 

•  Steps  to  solve  a  problem: 


•  Currently  requires  writing  a  C++  program 

•  Versatile  but  complex  approach 

•  Working  on  simpler  interface  for  SAMS 


Solve 

^  system  of 
equations 


DISTRIBUTION  STATEMENT  A:  Approved  for  public  release.  Distribution  is  unlimited. 


22 


Integration  with  SAMS 


•  Time-accurate  fluid-structure  interaction 

•  Multidisciplinary  adjoint  sensitivities 

•  Streamlined  problem  definition  and  execution 

•  Interface  with  geometry  definition 
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YAMSS 


Yet  Another  Mode-based  Structural  Solver 

•  Alternative  to  FUN3D  modal  solver 

•  Use  to  test  modularity  of  structural  solver 
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Test  Cases 
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Overview 


Case 

nD 

Flow 

Structure 

Vortex-induced  vibrations 

2D 

transient 

spring 

Airfoil  with  gust 

2D 

transient 

linear  beams 

Turek-Hron 

2D 

transient 

hyperelastic  flag 

AGARD  445.6 

3D 

steady 

linear  modes 

G-Wing 

3D 

transient 

aluminum  plate 

X-56A 

3D 

steady 

built-up  structure 

Tailless  supersonic  aircraft 

3D 

transient 

rigid 

Incrementally  demonstrate  capability 
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Vortex-Induced  Vibrations 


•  Cylinder  in  cross  flow,  free  to  translate  or  rotate 

•  Vortical  flow  may  induce  vibrations  in  cylinder 

•  Parameters:  k,  Re 


References 

Anagnostopoulos,  R,  and  P.  W.  Bearman.  “Response  characteristics  of  a  vortex-excited 
cylinder  at  low  Reynolds  numbers.”  Journal  of  Fluids  and  Structures  6  (1992):  39 — 50. 

Williamson,  C.  H.  K.,  and  R.  Govordhan,  “Vortex-induced  \/\t>ra\\or\s”  Annual  Review  of  Fluid 
Mechanics  36  (2004):  413—455. 
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Airfoil  in  a  Gust 


•  Sting-mounted  airfoil  in  a  transient  gust 

•  Parameters:  Ug^g^,  a,  p,  E 

I 


Reference 

Liu,  S.,  and  R.  Canfield,  “Continuum  shape  sensitivity  for  nonlinear  transient  aeroelastic  gust 
response.”  A\AA-20^^-^97^,  52^^ AIAA/ASME/ASCE/AHS/ASC  Structures,  Structural 
Dynamics  and  Materials  Conference,  2011. 

Cross,  D.  M.,  Local  continuum  sensitivity  method  for  shape  design  derivatives  using  spatial 
gradient  reconstruction.  Diss.  Virginia  Polytechnic  Institute  and  State  University,  2014. 
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Turek-Hron 


•  Rigid  cylinder  with  an  attached  elastic  structure 

•  Flow  induces  deformation  in  elastic  structure 

•  Parameters:  Re,  U,  p,  E,  p 


Reference 

Turek,  S.,  and  J.  Hron,  Proposai  for  numericai  benchmarking  of  fiuid-structure  interaction 
between  an  elastic  object  and  laminar  incompressibie  fiow.”  Fluid-structure  interaction. 
Springer  Beriin  Heidelberg,  2006,  371 — 385. 
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AGARD  445.6  Wing 


•  Standard  aeroelastic  benchmark  problem 


•  Modal  structure 

•  Parameters: p„ 


Reference 


< 


l9 


u 

h 

IL 

Eisoh 

HJ 

c 
>- 
G 

.-IQOh 


50  - 


tup. 

FUN3C-E 

FUN3D  NS  SA  QrkI 

FUM3D  HS  SA  Fine  Qrld 


'  NS  -NdvIh  Stakes 
^  Spalnrt-AJImara^. 


J_ 

0.4  O.S 


_L 

0.6 


I 


I 


o.T  o.a  0.& 

Macli 


1.1 


1.2 


Silva  et  a\.,AIAA,  2014-0496 


Yates,  Jr,  E.  C.,  “AGARD  standard  aeroelastic  configurations  for  dynamic  response. 

Candidate  configuration  I. -Wing  445.6.”  NASATM  100492  (1987). 

Lee-Rausch,  E.  M.,  and  J.  T.  Batina,  “Calculation  of  AGARD  wing  445.6  flutter  using  Navier- 
Stokes  aerodynamics.”  AIAA-1 993-3476,  1V^  Applied  Aerodynamics  Conference  (1993): 
9—11. 
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G-Wing 


•  Swept  transport  wing 


•  Nonlinear  structure,  exhibits  high-altitude  LCO 


•  Parameters:  M^,  h,  E 


Reference 


Bendiksen,  O.  O.  “Transonic  limit  cycle  flutter  of  high-aspect-ratio  swept  wings.”  Journal  of 
Aircraft  45.5  (2008):  1522—1533. 

Bendiksen,  O.  O.  “High-altitude  limit  cycle  flutter  of  transonic  wings.”  Journal  of  Aircraft  46. ^ 
(2009):  123—136. 
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X-56A 


•  Flight-tested  UAV  with  rigid  and  flexible  wings 

•  Flutter  suppression  &  gust  load  alleviation 


•  Parameters:  M^,  K 


Reference 


Beranek,  J.,  and  others.  “Conceptual  design  of  a  multi-utility  aeroelastic  demonstrator.”  AIAA- 
2010-9350.  13^^  A! AA/ISSMO  Multidisciplinary  Analysis  Optimization  Conference,  2010. 

Jones,  J.,  and  C.  E.  Cesnik.  “Nonlinear  aeroelastic  analysis  of  the  X-56  multi-utility 
demonstrator.”  AIAA-201 6-1 799.  15^'^  Dynamics  Specialists  Conference,  2016. 


DISTRIBUTION  STATEMENT  A:  Approved  for  public  release.  Distribution  is  unlimited. 


32 


Tailless  Aircraft 


•  Tailless  supersonic  aircraft 

•  Fuselage-mounted  tabs  improve  yaw  control 

•  Parameters:  tab  shape  and  position 


Meckstroth,  C.,  and  others.  “Aerodynamic  modeling  techniques  for  efficient  supersonic  air 
vehicle  multidisciplinary  design  optimization.”  AIAA-201 4-3254.  32'^^  AiAA  Applied 
Aerodynamics  Conference,  2014. 

Meckstroth,  C.,  and  W.  Blake.  “Control  focused  multidisciplinary  design  optimization  of  tailless 
fighter  aircraft.”  AIAA-201 5-2324.  1 6^*^  AIAA/ISSMO  Multidisciplinary  Analysis  and 
Optimization''  ' 
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Concluding  Remarks 
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Progress 


•  Major  changes  to  FUN3D  are  underway 

•  MAST  improvements  are  readying  it  for  SAMS 

•  Integration  infrastructure  ready  for  use 

•  Adopting  FUNtoFEM  (  for  FSI 


Next  step:  mode-based  static  aeroelasticity 
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Personnel 


•  SAMS  has  been  understaffed 

•  AFRL  has  had  two  unfilled  contractor  positions 

1 .  Filled  in  November 

2.  Currently  reviewing  applicants 

•  NASA  had  one  unfilled  position 

1 .  Filled  in  October 


Behind  schedule,  but  pace  accelerating 
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Digital  Collaboration 


•  Need  to  exchange  source  code  and  data  files 

•  Existing  ’’solutions”  not  working 

-  AFRL  unsuccessful  getting  full  access  to  LaRC 
GitLab 

-  NASA  unsuccessful  getting  accounts  to  AFRL  ELSZ 

•  Large  file  transfer  mechanisms  do  work 

-  But  not  good  for  collaboration 


Digitcd  collaboration  remains  a  significant  problem 
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